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The ferrofluid dynamics theory is applied to thermodiffusive problems in magnetic fluids in the 
presence of magnetic fields. The analytical form for the magnetic part of the chemical potential and 
the most general expression of the mass flux are given. By employing these results to experiments, 
global Soret coefficients in agreement with measurements are determined. Also an estimate for a 
hitherto unknown transport coefficient is made. 

PACS numbers: 66.10.Cb, 75.50.Mm, 47.10.+g 



I. INTRODUCTION 



Magnetic fluids (MFs) are colloidal suspensions of ferromagnetic nanoparticles dispersed in a nonmagnetic carrier 
liquid. MFs behave superparamagnetically in a magnetic field and have a far reaching application potential spanning 
from sealants in rotary shaft to heat dissipatcrs in loud speaker coils to carrier liquids for medical substances 
Starting in the mid-sixties of the last century, when MFs were first available, research on those fluids had been 
proceeding on the calm fairway of an established and well founded field of research. Particularly the theoretical work 
had been based on the achievements of the two pioneers, Rosensweig Q and Shliomis But fO years ago a series 

of papers 0, S B 13 started to appear pointing specifically to the deficiencies of the microscopic approach in '3, ^ and 
proposing "... a general, strictly macroscopic approach relying solely on symmetry considerations, conservation laws, 
and thermodynamics." This approach, called ferrofluid dynamics (FFD), sparked an impassioned discussion 

[Tol Ell about which theory explains better the experimental facts for the reduced viscosity of a MF in an ac magnetic 
field |f 2j or for the magnetovortical resonance (1^ J^J . That new theory also triggered an experiment confirming 
a proposed nonzero transport coefficient which is zero in the microscopic approach 0, . Other proposed effects as 
shear-excited sound p^ll7l| await their confirmation yet. 

For the description of thermal convection in magnetic fluid, the fluid has been considered as a one-component fluid 
with effective properties in many studies (see [18|, [lOj, (2^ and references therein) . The limits of this coarse grained 
view onto the colloidal suspension of ferromagnetic nanoparticles are just being revealed. Considering a magnetic 
fluid as a binary liquid, the thermal convection is found to set in at Rayleigh numbers well below the threshold for a 
MF considered as a single-component fluid ^2^]. 

The thermodiffusive or Soret effect describes the establishment of concentration gradients in response to temperature 
gradients for a two-(or multi-) component fluid. Since the motion of the ferromagnetic nanoparticles in the MF can be 
influenced by external magnetic fields, the Soret effect in MFs shows a strong dependence on any nonzero magnetic 
field strength 22^ 23, 2^ . In a vertical layer the Soret coefficient St depends raonmonotonously on the strength of 
the field in the cases where the field is either parallel or perpendicular to the temperature gradient [23 . |23| . Contrary, 
for both orientations of the magnetic field the Soret coefficient depends monotonously on the strength of the field if 
the layer is horizontal 0|. The changes of Sj can be up to six times its zero field value and even a change of 
the sign of St was measured for strong fields [22, 123, |2J| . 

The known theoretical approaches for the Soret effect in magnetic fluids [25l |2^ need as an essential input an 
expression for the magnetophoretic velocity of the nanoparticles with respect to the carrier liquid. For that purpose 
certain microscopic properties are assumed as a dilute colloid containing spherical particles of equal size and the 
applicability of the Stokes hydrodynamic drag psi |2^ . Also assumptions about the deformation of the temperature 
distribution around the particle are made if its thermal conductivity is different from that of the surrounding carrier 
liquid |2^. A comparison with the known experimental results shows great differences: the microscopic theory pMI?^ 
gives only changes of St which are about three orders of magnitude smaller than the experimentally measured ones 
([23II and Fig. 23 in [l^). Also in the frame of a thermodynamic approach it is not possible to describe the drastic 
changes of St measured in the experiment. That means that with respect to thermodiffusive processes in MF in the 
presence of magnetic fields a wide gap between experiment and theory has to be bridged. Therefore it is the aim of 
this work to present a different approach, where in the frame of a macroscopic theory, the FFD, the experimental 
results can be described significantly better. 

Usually an external temperature gradient causes both convection and thermodiffusion in any colloidal suspension. 
How these two effects are interacting with each other is not yet finally resolved as the discussion about the possibility 
of a state of relaxation-oscillation convection highlights |2l|, [23 . The mutual interference of convection and thermodif- 



fusion is even more severe if additionally an external magnetic field is applied as in the case of MFs |22ll23ll24ll26l l3 
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The problems caused by that mutual interference for the determination of the Soret effect are outlined in |2J| and 
result in a new experimental setup for a horizontal layer of MF which is analyzed theoretically in this work. 

II. FERROFLUID DYNAMICS: CHEMICAL POTENTIAL AND MASS FLUX 

The macroscopic FFD approach is presented without magnetodissipation, i.e. the magnetization M is always 
parallel to the magnetic field H, but with dissipative mass fluxes for the two constituents of the MF. The analysis 
will result in an analytical expression for the magnetic part of the chemical potential and a general expression for the 
mass flux without any assumption about the properties of the MF and the temperature distribution. 

The principal structure of the ferrofluid dynamics theory was given in j£j . It is based firstly on general principles 
as symmetry considerations and conservations laws and on irreversible thermodynamics. The second independent 
component of which a macroscopic theory is made of is the set of material-dependent parameters as susceptibilities 
and transport coefficients. The latter can be determined by suitable experiments which are used here to determine 
transport coefficients for thermodiffusive processes in magnetic fluids in the presence of magnetic fields. 

As usual in theories based on thermodynamical considerations, one starts with the thermodynamic energy density 
u. It is taken as a function of the entropy density s, the density p*^^^ of the magnetic part of the fluid, the momentum 
density g = pv, the total density p, and the magnetic induction B = /io(M + H) 0, 

du = Tds + jlcdp^'^^ + Vidgi + ^x^'^^dp + H,dBi , (2.1) 

where pc = A^^'' ~ M*"^-* is the difference in the chemical potentials of the two constituents. The conservation laws for 
the density of the magnetic and nonmagnetic part p*-^) are 

^ -V.(p(%.-jf) , (2.2) 
= -V,(p(^)..+jf) , (2.3) 

where jf^^^ — —ji""^^ —jf was used to ensure the conservation of the total density p = p^^^ + /9^^^=</)pm + (1 ^ <f>)Pc\- 
The density of the magnetic particles (carrier liquid) is denoted by Pm (Pci) and 4> is the volume fraction of magnetic 
particles in the fluid. The dissipative mass flux is proportional to the gradient of the chemical potential with 
Ac = P'cip, P^^^KTjV,!!) and the temperature gradient j^. It is assumed that the magnetic part of the chemical 
potential can be separated (s^. 

p, = Pcip, p(i) , T, v) + pr (P, P^'^ , T, V, H) . (2.4) 

This assumption guarantees a nonzero chemical potential for H = and is confirmed by calculations for MF with 
chains, where the magnetic part contributes additive to the total chemical potential |33j . The nonmagnetic part of 
the chemical potential is given by pc — (kBT/m^) Inci — (k^T / mc\) lnc2, where mm (wd) is the mass of a magnetic 
(carrier liquid) particle Isj- 

The experiments 123 . |23| show that any nonzero strength of the magnetic field influences the thermodiffusive 
processes. Thus the general ansatz for the dissipative mass flux is (following the notation in Q) 

]? = 6V,r + Sy.pc + Cii A/,Mj Vj/ic + Cxe^jfcMj VfePe . (2.5) 

Whereas the flrst two terms characterize isotropic mass fluxes caused by gradients in the temperature or in the chemical 
potential, the last two terms describe anisotropic mass fluxes, namely parallel and perpendicular to the direction of 
M. The last term corresponds to that one in the analogous ansatz for the heat flux, where the phenomenon is called 
transversal Righi-Leduc effect |3lj| since the primary current is perpendicular to the produced effect. 

It was emphasized in that the "proliferation of transport coefficients" , i.e. ^ — > (^, , ), takes place in the case 
of strong magnetic fields. But the experiments show that at least for thermodiffusive processes that general statement 
seems not to be true. In the figures in [23, IH, 01 with respect to the changes of St it is evident that small magnetic 
fields in the order of less than 50 kA/m are sufficient to generate effects, where one can clearly distinguish between a 
parallel or a perpendicular orientation between temperature gradient and field. Therefore the in ^ firstly introduced 
coefficients ^|| and are considered here as nonzero for all magnetic field strengths. 

With the above given dependences of the chemical potential in Eq. H2.4I) , its gradient is 



(2.6) 
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The first expression in Eq. (|2.6|) will become later the term for the barodifFusion and can be neglected in an incom- 
pressible fluid not subjected to any pressure gradient. For the fourth and fifth term [28Ll3^ hold 



The transformation u = u — vjgj — HjBj was made in order to match the dependences of the energy density and the 
chemical potential and usage of the fact that derivatives of quantities are zero which are independent of each other. 
From the last equality the analytical result for the magnetic part of the chemical potential follows, 

DM 

Z^™--^"/ ^rfi?', (2.9) 

where M and H denote the absolute value of the magnetic field and the magnetization. Eq. (|2.9|) allows a direct 
determination of if the magnetization M{H, p*-^-* , T) is known without any assumption about the properties of the 
MF in contrast to |27l l2a. 133 . l34j . According to these references the determination of the chemical potential needs the 
knowledge of quantities like the volume concentration of the nanoparticles ,27. .pS.! or the strength of the magnetodipole 
interaction [33 or the effective field experienced by a single particle in the MF |23,|2a|33|- Compared with the effort to 
evaluate these microscopic details, the advantage of the macroscopic approach of the FED is apparent. A measurement 
of the magnetization as function of the magnetic field and the density is sufficient to determine the chemical potential 
for any magnetic fluid. 

Inserting Eq. 1)2. 6|l into Eq. H2.5|l and using Eqs. (|2.7M2.8f) an expression for the mass flux results. 



p p dT J p dT dT 



-^M (MVT) + — (M X VT) 
P P 



-^^[eVi/ + e||M(MVif)+ex(MxVi/)] , (2.10) 

which is generally valid, independent of the size distribution of the magnetic particles, concentration inhomogeneities 
in the suspension or the form of the temperature gradient. Therefore Eq. (|2.10|l is the generalization of the mass flux 
given in |25l |. The concentration of the magnetic particles ci — p^^^ j p is defined by means of the mass fraction of 
the total density p (3ll|. The first four terms describe mass flow caused by thermophoresis (~ VT), the second four 
terms by diffusiophoresis (~Vci), and the last three by magnetophoresis (~ Vi?). There are two unknown transport 
coefficients, ^|| and ^x, since for zero magnetic field, Eq. (|2.10() reduces to the classical result (see Eq. (227), Chapt. 
XI in 31]) 

(i\ , i dp.c\ , . dp., 



P ' pdTj^^^^W^^^^ 
= CiC2DtVT + D^ci (2.11) 

with {Dt) Dc the (thermal) diffusion coefficient known for MFs from previous experiments Js^ and C2 = 1 — ci. 
According to the philosophy of the FFD approach, the determination of the unknown transport coefficients ^|| and 
needs suitable experiments which were conducted just recently (24| . 



III. APPLICATION TO EXPERIMENTS AND DISCUSSION 

According to the experiments for a horizontal layer of MF of thickness h , a horizontally unbounded layer of a 
dielectric, viscous, and incompressible MF sandwiched between two perfect conducting plates is considered. The lower 
plate is cooled to Ti and the upper one is heated to T2. The resulting temperature gradient stabilizes the quiescent 
conductive state. From the equation of heat conduction. 



(3.1) 
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and the boundary conditions, 

T{z - h/2) - T2 
the temperature profile of the conductive state 



T[z = -h/2) = Ti 



(3.2) 



(3.3) 



follows with T{) = (Ti + T2)/2 and k denotes the thermal diffusivity. Since the plates are impenetrable, the diffusion 
equation, 



dci (} 



has to be supplemented with the boundary condition 

jf(z = ±V2) = 0. 
Rearranging this boundary condition with the help of Eq. (|2.11|) , 



-h 



dci 



ciC2(T2 - Ti) dz 



z=±h/2 



Dt 



(3.4) 



(3.5) 



(3.6) 



the Soret coefhcient in the zero field case can be expressed. In the same way the global Soret coefficient, measured in 
[2^. in the presence of a magnetic field can be determined by using Eq. (|2.10() . 

If a spatially homogeneous static magnetic field is applied perpendicular to the layer, the resulting magnetic field 
gradient inside the fluid is parallel to the temperature gradient. Therefore this setup is called parallel and is analyzed 
first. 

Taking diffusion processes into account, the magnetization in the fluid can be written in the form 



M = 



Mo + x{H - Ho) - K{T - To) 



dM 



- M 



(3.7) 



where Mq = Mq^HqjTq, (po) is the reference magnetization belonging to the reference values Hn, Tq, and (jjo for the 
magnetic field, the temperature, and the volume fraction. Extending the expressions given in [l^, magnetization and 
magnetic field for the conductive state are 



M 

H = Ho- 



hil + x) 
K{Ti - T2) 

h{l + x) 



z - N{ci - ci,o) 



(3.8) 
(3.9) 



with the susceptibility x = dM/dH, the pyromagnetic coefficient K = —dM/dT, the dcnsomagnetic coefficient 
N = dM/dci = [p/ prn){dM/d4>\ and cis) = ci(To). Inserting Eqs. (EH hito Eq. (|TTUI) and rearranging the 
boundary condition l|3.5() in the same manner as in the zero field case, the global Soret coefficient in the parallel setup 
reads: 



St + ^ 

I _ Cl.»C2fiP 



Dc dT 



Dc dT 



Po 



dM 



Dc Dc y(i + x) 



e dp^ , ^11 dpc 



Dc dpW Dc dpW 



po dM 
'J dpW 



1. + ^Af2 
Dc Dc 



(3.10) 



N 



Knowing M{H, p^^\T) in analytical form allows one to calculate p^ and its derivatives. The measured magnetiza- 
tion curve (from Fig. 57 in f22j) could be nicely fitted with M = MhX^(pL{Xda), where L(Ada) = coth(Ada) — l/{XdC() 
is the Langevin function, a = po'mHl[k^T) the Langevin parameter, m — MbTr(P/6 the magnetic moment of a 
particle, and the Boltzmann constant. A^; and A^ are two geometrical fit parameters. They refiect small deviations 
from the volume fraction (p = 0.2 and the (5-shaped size distribution (Fig. 59 in [l^). Using A^ = 0.99, A^ = 0.84 and 
the material data Mb = 450 kA/m (magnetization of the magnetic bulk solid), d = 9 nm, pm = 5.15 g/cm'^ from [2^ . 
the solid line in Fig. ^ shows a very good agreement with the measured magnetization (o). Considering the chosen 
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values for and A^, only the volume fraction cj) had to be adjusted to the measured data. Variations in (j) are likely 
caused by a nonmagnetic surface layer of the nanoparticles |3.6| and its solvability in the carrier liquid. According to 
the statement at the beginning of the paragraph, one has 



{In [coth^Xaa) - l] + 21n(Ada)} 



Ad 2pmTO 

from which one can calculate the derivatives with respect to T and p^^^ 



(3.11) 




100 200 300 

magnetic field H (kA/m) 



FIG. 1: Experimental data (o) from Fig. 57 in |23l and theory (solid line) for the magnetization at room temperature T = 293 
K. The details of the used Langevin function are given in the text. 



With the pyromagnetic coefRcient K taken from Fig. 4 in J22|, it remains the four unknowns St, Dc, and ^|| 
in Eq. H3.10() . To fit s!^ to the experiment, the combined quantities ^/Dc and ^\\/Dc arc used as fit parameters, 
since 5*7- = 0.15 was measured in the zero field case 24] but not Dc- The solid line in Fig. [21gives the best two 
parameter fit, yielding ^/I?c = 8.2 kgs^/m^ and ^\\/Dc = —lAl ■ 10^'' kgs^/(m'^ A^). The difference in the absolute 
values of about 8 orders of magnitude is not surprisingly, since one would assume such a relation according to the 
argument that anisotropic fluxes in the mass flux (|2.5|l are relevant only for strong fields 9]. Inspecting Eq. (|3.10() 
closer, it is revealed that ^\\/Dc is multiplied by which gives already for small magnetic fields a factor of ^^10^. 
The two other terms are not so relevant because < d^c/dT < 0.016 J/(K kg) and < dfic/dp'^^^ < 0.35 Jm^/kg^ 
for < H < 350 kA/m. To underline the relevance of £,\\/Dc even for small fields, the dot-dashed line in Fig. [21 

displays s!^ for ^\\/Dc — and all other parameters as before. Now the theoretical curve misses the measured data (□) 
clearly. Taking a typical value for the diffusion coefficientj_Dc ~ 10^^^ m^/s the new transport coefficient can be 
estimated to ^|| ~ — 10~^* kgs/(mA^) for the MF in j22.l24l |. Thus those experiments deliver the necessary input for 
determining the material-dependent transport coefficients which are a priori unknown in a macroscopic theory as the 
FFD. Another example for the experimental determination of diffusion and thermodiffusion coefficients is presented 
in (35j , whereas in 37] these coefficients were calculated on the basis of a microscopic theory. 

In contrast to the parallel setup, in the perpendicular setup the spatially homogeneous static magnetic field is 
applied perpendicular to the temperature gradient, i.e. the magnetic field is parallel to the layer. The diffusion 
equation gets now the form 



dci 



dp. 



H 

Qy2 Q^2 



(3.12) 



The boundary condition for the z-component of the mass flux yields 



dci 
dz 



dci 
? dy 



Ci C2 Dt 



p dT 



at z = ±h/2 . 



(3.13) 
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FIG. 2: Global Soret coefficients Sj, and 5*^ against the magnetic field strength for the parallel {H \\ VT) and perpendicular 
setup {HIVT). The solid line shows the best fit of [see Eq. (ICTil l with ^/D^ = 8.2 kgs^/m'' and C||/-Dc = -1.41 ■ 10"'' 
kgs^/(m'' A^) to the experimental data (□). The dot-dashed line displays S^^ for the same parameters but £,\\/Dc = 0. The 



dashed line indicates the best fit of [see Eq. lETT^ I with ^/D^ = 8.2 kgsVm^ and F = 3.75 ■ 10" 
experimental data (A). For F and all other values see text. 



kgs^/(m^A) to the 



Since no analytical solution for that boundary value problem is known, the following coarse approximation is made: 
{dci/dy)z=±h/2 shall be a constant C for all iJ-values tested here. The global Soret coefficient in the perpendicular 
setup can be approximated than by 



St 



e 1 



where F = {^±C)/Dc will be used as the only fit parameter since £,/Dc was determined in the parallel setup. With 
-Ti = l K, = 0.2, and h^l mm (2^, the best fit yields i^ = 3.75 • 10"^ kgs2/(m5A). The inferior match with the 
experimental data (see A and dashed line in Fig. I^J in comparison with the parallel setup is due to the approximation 
that idci/dy)z=±h/2 is constant. In the real system it will depend on the magnetic field since the solution for ci 
depends on the magnetic field. 



IV. CONCLUSION 



The ferrofluid dynamics theory is applied to thermodiffusive problems in magnetic fluids in the presence of magnetic 
fields, where the MF is considered as a binary mixture. In the frame work of this theory the chemical potential could 
be determined analytically. Also a general expression for the mass flux is given which is independent of the fluid 
properties, temperature distribution and assumptions about the concentration of the nanoparticles. Applying these 
results to the experiments |23|, their data could be interpreted better (see Fig.|2l than with the previous theory 
which gave values about three orders of magnitude too small. Three transport coefficients, which are inherent parts the 
macroscopic ferrofluid dynamics theory 9] , had to be used to fit this theory with the only sets of experiments available 
at present. In general, it is shown that for thermodiffusive problems in magnetic fluids, i.e. in colloidal suspensions 
sensitive to external fields, anisotropic mass fluxes are relevant and no small contributions for any nonzero strengths 
of the magnetic field. To elucidate this insight, more well designed experiments and further theoretical analyses are 
needed to improve the knowledge about thermodiffusive processes in magnetic fluids. 
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